Clustering, Angular Size and Dark Energy 
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The influence of dark matter inhomogeneities on the angular size-redshift test is investigated for 
a large class of flat cosmological models driven by dark energy plus a cold dark matter component 
(XCDM). The results are presented in two steps. First, the mass inhomogeneities are modeled 
by a generalized Zeldovich-Kantowski-Dyer-Roeder (ZKDR) distance which is characterized by a 
smoothness parameter a(z) and a power index 7, and, second, we provide a statistical analysis to 
angular size data for a large sample of milliarcsecond compact radio sources. As a general result, we 
have found that the a parameter is totally unconstrained by this sample of angular diameter data. 

PACS numbers: 98.80.-k; 95.36.+x; 95.35.+d 



I. INTRODUCTION 



An impressive convergence of recent astronomical 
observations are suggesting that our world behaves like 
a spatially flat scenario dominated by cold dark matter 
(CDM) plus an exotic component endowed with large 
negative pressure, usually named dark energy [j], 0, Ell- 
in the framework of general relativity, besides the cos- 
mological constant, there are several candidates for dark 
energy, among them: a vacuum decaying energy density, 
or a time varying AJt) Q, the so-called "X-matter" Q, 
a relic scalar field @, and a Chaplygin Gas Some 
recent review articles discussing the history, interpreta- 
tions, as well as, the major difficulties of such candidates 
have also been published in the last few years Q. 



In the case of X-matter, for instance, the dark energy 
component is simply described by an equation of state 
p x = Lop x . The case uj = — 1 reduces to the cosmo- 
logical constant, and together the CDM defines the sce- 
nario usually referred to as "cosmic concordance model" 
(ACDM). The imposition u> > —1 is physically motivated 
by the classical fluid description [fj]. However, as dis- 
cussed by several authors, such an imposition introduces 
a strong bias in the parameter determination from obser- 
vational data. In order to take into account this difficulty, 
superquintessence or phantom dark energy cosmologies 
have been recently considered where such a condition is 
relaxed [n| ■ In contrast to the usual quintessence model, 
a decoupled phantom component presents an anomalous 
evolutionary behavior. For instance, the existence of fu- 
ture curvature singularities, a growth of the energy den- 
sity with the expansion, or even the possibility of a rip-off 
of the structure of matter at all scales are theoretically ex- 
pected ([ll[ for a thermodynamic discussion). Although 
possessing such strange features, the phantom behavior 
is theoretically allowed by some kinetically scalar field 
driven cosmology as well as, by brane world mod- 
els [l3j], and, perhaps, more important to the present 
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work, a PhantomCDM cosmology provides a better fit to 
type la Supernovae observations than does the ACDM 
model 14]. Many others observational and theoretical 
properties phantom driven cosmologies (more generally, 
of XCDM scenarios) have been successfully confronted to 
standard results (see, for instance pH [Tr3. [l7l [l8l ) . 

In this context, one of the most important tasks for 
cosmologists nowadays is to confront different cosmolog- 
ical scenarios driven by cold dark matter (CDM) plus a 
given dark energy candidate with the available observa- 
tional data. As widely known, a key quantity for some 
cosmological tests is the angular distance-redshift rela- 
tion, Da{z), which for a homogeneous and isotropic back- 
ground, can readily be derived by using the Einstein field 
equations for the Friedmann-Robertson- Walker (FRW) 
geometry. From Da(z) one obtains the expression for 
the angular diameter 9{z) which can be compared with 
the available data for different samples of astronomical 
objects [20l | . 

Nevertheless, the real Universe is not perfectly homo- 
geneous, with light beams experiencing mass inhomo- 
geneities along their way. Actually, from small to in- 
termediate scales (< lOOMpc), there is a lot of structure 
in form of voids, clumps and clusters which is probed 
by the propagating light [2l[. Since the perturbed met- 
ric is unknown, an interesting possibility to account for 
such an effect is to introduce the smoothness parameter a 
which is a phenomenological representation of the mag- 
nification effects experienced by the light beam. From 
general grounds, one expects a redshift dependence of a 
since the degree of smoothness for the pressureless mat- 
ter is supposed to be a time varying quantity [13, OH- 
When a = 1 (filled beam), the homogeneous FRW case 
is fully recovered; a < 1 stands for a defocusing effect 
while a = represents a totally clumped universe (empty 
beam). The distance relation that takes these mass in- 
homogeneities into account was discussed by Zeldovich 
[22| followed by Kantowski [23| , although a clear-cut ap- 
plication for cosmology was given only in 1972 by Dyer 
and Roeder (24[. Later on, by considering a perturbed 
Friedmannian model Tomita [29j performed N-body sim- 
ulations with the CDM spectrum in order to determine 
the distribution for a (see also Ref. [3(| for a more gen- 
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eral analysis involving distances in perturbed models). 
Many references may also be found in the textbook by 
Schneider, Ehlers and Falco, [2^|, as well as, in Kantowski 

Many studies involving the ZKDR distances in dark en- 
ergy models have been published in the literature. Ana- 
lytical expressions for a general background in the empty 
beam approximation (a — 0) were derived by Sereno et 
al. [3l|. By assuming that both dominant components 
may be clustered they also discussed how the critical 
redshift, i.e., the value of z for which Da(z) is a max- 
imum (or 0(z) minimum), and compared to the homo- 
geneous background results as given by Lima and Alcaniz 
|32| , and, further discussed by Lewis and Ibata [Hj] , and 
Arau.jo and Stoeger [34[. More recently, Demianski et 
al. [35|, derived an useful analytical approximate solu- 
tion for a clumped concordance model (ACDM) valid on 
the interval < z < 10. Additional studies on this sub- 
ject involving time delay (Lewis and Ibata [33j] ; Giovi 
and Amendola [36|), gravitational lensing (Kochanek; 
Kochanek and Schechter [37]]) or even accelerated mod- 
els driven by particle creation have also been considered 

Although carefully investigated in many of their the- 
oretical and observational aspects, an overview in the 
literature shows that a quantitative analysis on the influ- 
ence of dark energy in connection with inhomogeneities 
present in the observed universe still remains to be stud- 
ied. Analytical expression for a general applied for the 
9(z) statistics with basis on a ACDM cosmology with 
constant a 40] . It was concluded that the best fit model 
occurs at Cl^j = 0.2 and a = 0.8 whether the characteris- 
tic angular size I of the compact radio sources is marginal- 
ized. More recently, the smoothness a parameter was 
constrained through a statistical analysis involving Su- 
pernovae la data [4l|. A ^-analysis based on the 182 
SNe la data of Riess et al. [2J constrained the pair of 
parameters (flM,ce) to be Qm = 0.33l^o7 an( ^ a — 0-42 
(2er). Such an analysis has also been carried out in the 
framework of a ACDM cosmology. 

In this paper, we focus our attention on X-matter 
cosmologies with special emphasis to phantom models 
(lj < —1) by taking into account the presence of a clus- 
tered cold dark matter. The mass inhomogeneities will 
be described by the ZKDR distance characterized by a 
smoothness parameter a(z) which depends on a positive 
power index 7. The main objective is to provide a sta- 
tistical analysis to angular size data from a large sample 
of milliarcsecond compact radio sources [42|] distributed 
over a wide range of redshifts (0.011 < z < 4.72) whose 
distance is defined by the ZKDR equation. As an extra 
bonus, it will be shown that a pure CDM model (Slj\/ = 1) 
is not compatible with these data even for the empty 
beam approximation (a = 0). 

The manuscript is organized as follows. In section 2 
we outline the derivation of the ZKDR equation for a X- 
CDM cosmology. We also provide some arguments (see 
Appendix) for a locally nonhomogeneous Universe where 



the homogeneous contribution of the dark matter obeys 
the relation ph — ap (pM I Po) 1 where 7 is a positive num- 
ber, pm is the average matter density and p its present 
value. In section 3 we analyze the constraints on the free 
parameters a and f2j\/ from angular size data. We end 
the paper by summarizing the main results in section 4. 



II. THE EXTENDED ZKDR EQUATION 

Let us now consider a flat FRW geometry (c = 1) 

ds 2 = dt 2 - R 2 {t)(dr 2 +r 2 d0 2 + r 2 sin 2 9 d(/) 2 ) , (1) 

where R(t) is the scale factor. Such a spacetime is sup- 
ported by the pressureless CDM fluid plus a X-matter 
component of densities pu and p x , respectively. Hence, 
the total energy momentum tensor, T^ v = T^ v ( M ) + 
T^ u m, can be written as 



T"" = [p M + (1 + u) Px ]U^U v - u>p x g^ } 



(2) 



where = 6% is the hydrodynamics 4-velocity of the 
comoving volume elements. In this framework, the in- 
dependent components of the Einstein Field Equations 
(EFE) 
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take the following forms 
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- = — Hf, 
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n M (^) 3 + (3^ + i)^(^) 3 ^ 



(3) 



(4) 



(5) 



where an overdot denotes derivative with respect to time 
and H a = 100/iKms _1 Mpc _1 is the Hubble parameter. 
By the flat condition, Q x = 1 — J1m, is the present day 
dark energy density parameter. As one may check from 
(2)-(5), the case u> = —1 describes effectively the favored 
"cosmic concordance model" (ACDM). 

On the other hand, in the framework of a comformally 
flat FRW metric, the optical scalar equation in the ge- 
ometric optics approximation reads (Optical shear ne- 
glected) [43] 



(6) 



where A is the beam cross sectional area, plicas means 
derivative with respect to the affine parameter describ- 
ing the null geodesies, and fc M is a 4- vector tangent to the 
photon trajectory whose divergence determines the opti- 
cal scalar expansion |TtI l3lL Ha] ■ The circular frequency 
of the light ray as seen by the observer with 4-velocity 
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U a is lu — U a k a , while the angular diameter distance, 
Da, is proportional to V^4 [25[. 

As widely known, there is no an acceptable averag- 
ing procedure for smoot hing out local inhomogeneities. 
After Dyer and Roeder [24(j . it is usual to introduce a 



phcnomcnological parameter, a(z) = 1 — 



Pel 

<Pm> 



, called 



the "smoothness" parameter. For each value of z, such 
a parameter quantifies the portion of matter in clumps 
(pel) relative to the amount of background matter which 
is uniformly distributed (pm). As a matter of fact, such 
authors examined only the case for constant a, however, 
the basic consequence of the structure formation process 
is that it must be a function of the redshift. Combining 
equations (2), (3) and (6), after a straightforward but 
lengthy algebra one finds that the angular diameter dis- 
tance, Da{z), obeys the following differential equation 



d 2 D A 



+ {1 + z)G—^- +HD A = 0, 



dz 2 ' " dz 

which satisfies the boundary conditions: 

D A (0) = 0, 



(ID A I _ 1 

~dT\o - l - 



(8) 



The functions T , Q and TL in equation (7) read 




Redshift (z) 



FIG. 1: The smoothness parameter as a function of the red- 
shift for some selected values of f3 and 7. All curves ap- 
proach the filled beam result (a = 1) at high redshifts re- 
gardless of the values of j3 and 7. Note that j3 determines 
a = a(z = 0). For a given j3 the curves starts at the same 
point but the rate approaching unit (filled beam) depends on 
the 7 parameter. 
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n M (i + z) 3 + 
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(9) 



The smoothness parameter a(z), appearing in the ex- 
pression of TL, assumes the form below (see Appendix A 
for a detailed discussion) 



a(z) — 



/3o(l + *) 37 
1 + /3 (1 + z) 3 t' 



(10) 



where (3 and 7 are constants. Note that the fraction 
c«o = Po/ (1+A>) is the present day value of a(z). In Fig. 
1 we show the general behavior of a(z) for some selected 
values of f3 and 7. 

At this point, it is interesting to compare Eq. (7) 
together the subsidiary definitions (8)-(10) with other 
treatments appearing in the literature. For 7 = (con- 
stant a) and w = — 1 (ACDM) it reduces to Eq. (2) as 
given by Alcaniz et al. [4(J. In fact, for u> = — 1 the 
function TL is given hy TL = ^J2a/(1 + z) 3 . Further, 
recalling the existence of a simple relation between the 
luminosity distance, and the angular-diameter distance 
(from Etherington principle [44], Dl = (1 + z) 2 Da), it is 
easy to see that Eq. (3) of Santos et al. [4l| is recovered. 



A more general expression for ACDM model (by includ- 
ing the curvature term) has been derived by Demianski et 
al. [35j | - As one may check, for a constant, by identifying 
u> = m/3 — 1, our Eq. (7) is exactly Eq.(10) as presented 
by Giovi and Amendola [36] in their time delay studies 
(see also Eq. (2) of Sereno et al. [45j]). Different from 
other approaches appearing in the literature (see for in- 
stance, Refs. [1^, Ho]), we stress that in this paper the a 
parameter is always smaller than unity. In addition, the 
a parameter may also depend on the direction along the 
line of sight (for a discussion of such effects see Linder 
[3, Sereno et al. [45], Wang H). 

Let us now discuss the integration of the ZKDR equa- 
tion with emphasis in the so-called phantom dark energy 
model (lo < — 1). In what follows, assuming that a; is a 
constant, we have applied for all graphics a simple Runge- 
Kutta scheme (see, for instance, the rksuite package from 
www.netlib.org) . 

In Figure 2 one can see how the equation of state pa- 
rameter, co, affects the angular diameter distance. For 
fixed values of Qm = 0.3, [3 a = 0.5 and 7 = 0, all the dis- 
tances increase with the redshift when uj diminishes and 
enters in the phantom regime (oj < —1). For compari- 
son we have also plotted the case for ACDM cosmology 
(a, = -l). 

In Fig. 3 we show the effect of the 7 parameter on the 
angular diameter distance for a specific phantom cosmol- 
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- in = -1.5 (Phantom) 
1.4 (Phantom) 
in = -1.3 (Phantom) 
!) = -1.2 (Phantom) 
I (Phantom) 
a = -1.0 (ACDM) 



i 0.5 i = o n = 0.3 
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FIG. 2: Angular diameter distance for a flat FRW phantom 
cosmology. The curves display the effect of the equation of 
state parameter for (5 — 0.5 and 7 = 0. The thick curve 
corresponds to the ACDM model. Note that for a given red- 
shift, the distances always increase for uj beyond the phantom 
divide line (to < —1). 



ogy with uj = — 1.3, as requested by some recent analyzes 
of Supernovae data 0] . For this plot we have considered 
(3 = 0.5. As shown in Appendix A, /3 Q = {ph/ Pci)z=o, is 
the present ratio between the homogeneous (ph) and the 
clumped (p c i) fractions. It was fixed in such a way that 
a assumes the value 0.33. Until redshifts of the order 
of 2, the distance grows for smaller values of 7, and after 
that, it decreases following nearly the same behavior. 

In Fig. 4 we display the influence of the [3 Q parameter 
on the angular diameter distance for two distinct sets of 
7 values. The cosmological framework is defined £Im = 
0.3 and the same equation of state parameter u> = —1.3 
(phantom cosmology). For each branch (a subset of 3 
curves with fixed 7) the distance increases for smaller 
values of (3 Q , as should be expected. 



III. ZKDR DISTANCE AND ANGULAR SIZE 
STATISTICS 

As we have seen, in order to apply the angular diame- 
ter distance to a more realistic description of the universe 
it is necessary to take into account local inhomogeneitics 
in the distribution of matter. Similarly, such a statement 
remains true for any cosmological test involving angular 
diameter distances, as for instance, measurements of an- 
gular size, 9(z), of distant objects. Thus, instead of the 
standard FRW homogeneous diameter distance one must 



FIG. 3: Effects of the 7 parameter on the angular diameter 
distance. For all curves we fixed u> = —1.3, /3 = 0.5 and 
= 0.3. Note that the distances increase for smaller values 

of 7. 



consider the solutions of the ZKDR equation. 

Here we are concerned with angular diameters of light 
sources described as rigid rods and not isophotal diam- 
eters. In the FRW metric, the angular size of a light 
source of proper length I (assumed free of evolutionary 
effects) and located at redshift z can be written as 



9{z) 



D A (z) 



(11) 



where I = lOOZh is the angular size scale expressed in 
milliarcsecond (mas) while I is measured in parsecs for 
compact radio sources (see below). 

Let us now discuss the constraints from angular size 
measurements of high z objects on the cosmological pa- 
rameters. The present analysis is based on the angu- 
lar size data for milliarcsecond compact radio sources 
compiled by Gurvits et al. [UJ (see also for ap- 
plications to the unclustered FRW case). This sample 
is composed by 145 sources at low and high redshifts 
(0.011 < z < 4.72) distributed into 12 bins with 12- 
13 sources per bin (for more details see Gurvits et al. 
[42l|). In Figure 5 we show the binned data of the me- 
dian angular size plotted as a function of redshift z to 
the case with 7 = and some selected values of VIm and 
a Q = /?o/(l — Po) — constant. As can be seen there, for 
a given value of Q.m the corresponding curve is slightly 
modified for different values of the smoothness parameter 
a. 

Now, in order to constrain the cosmic parameters, we 
first fix the central value of the Hubble parameter ob- 
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FIG. 4: Influence of the j3 parameter on the angular diameter 
distance for Qm = 0.3 and u) = —1.3. The curves are sepa- 
rated in two sets corresponding to the values of 7 = 0.5, 0.9 
as indicated in the box. As expected, both sets present the 
same behavior at low redshifts. 



FIG. 5: Angular size versus redshift according to the ZKDR 
distance. Curves for Qm = 0.3, 7 = and different values 
of lu are shown. The data points correspond to 145 compact 
radio sources binned into 12 bins (Gurvits et al. [42J). For 
comparison the filled beam ACDM has been included. 



tained by the Hubble Space Telescope (HST) key project 
H = 72±8km.s" 1 .Mpc" 1 (Freedman et al. [13). Nowa- 
days, this HST result has been confirmed by many dif- 
ferent classes of estimators like the Sunyaev-Zeldovich ef- 
fect and the ages of old high redshifts galaxies [48| . This 
value is also in accordance with the 3 years release of the 
WMAP team [3j, however, it is greater than the recent 
determination by Sandage and collaborators [4!| . Follow- 
ing standard lines, the confidence regions are constructed 
through a x 2 minimization 

2/, n [9(zi,l,uj,a) ~ m ] 2 
X (i.w.a) = 2^ 2 ' ( 12 ) 

where 0(zi, I, ui, a) is defined from Eq. (7) and O i are 
the observed values of the angular size with errors <Ji 
of the ith bin in the sample. The confidence regions 
are defined by the conventional two-parameters x 2 levels. 
In this analysis, the intrinsic length I, is considered a 
kind of "nuisance" parameter, and, as such, we have also 
marginalized over it. 

In Fig. 6 we show confidence regions in the to — a plane 
fixing Qm = 0.263, and assuming a Gaussian prior on the 
lo parameter, i.e., lo = — 1 ± 0.3 (in order to accelerate 
the universe). The "x" indicates the best fit model that 
occurs at ui = —1.03 and a ~ 0.9. 

In Fig. 7 the confidence regions are shown in the 
Qm — a plane. We have now assumed a Gaussian prior 
on Qm, i.e., Qm = 0.3 ±0.1 from the large scale struc- 
ture. From Figs. 6 and 7, it is also perceptible that 



while the parameters to and Qm are strongly restricted, 
the entire interval of a is still allowed. This shows the 
impossibility of tightly constraining the smoothness pa- 
rameter a with the current angular size data. This re- 
sult is in good agreement with the one found by Lima 
and Alcaniz 0] where the same data set were used to 
investigate constraints on quintessence scenarios in ho- 
mogeneous background, and is also in line with the one 
obtained by Barber et al. [5(| who argued in favor of 
a Q = a{z = 0) near unity (see also Alcaniz, Lima and 
Silva [4f| for constraints on a clustered ACDM model). 



IV. SUMMARY AND CONCLUDING 
REMARKS 

All cosmological distances must be notably modified 
whether the space-time is filled by a smooth dark energy 
component with negative pressure plus a clustered dark 
matter. Here we have addressed the question of how 
the angular diameter distance of extragalactic objects 
are modified by assuming a slightly inhomogeneous uni- 
verse. The present study complements our previous re- 
sults [2(3] by considering that the inhomogeneities can be 
described by the Zeldovich-Kantowski-Dyer-Roeder dis- 
tance (in this connection see also, Giovi and Amendola 
[3l| ; Lewis and Ibata [33| ; Sereno et al. [45| ; Demianski et 
al. (35}). The dark energy component was described by 
the equation of state p x = 10 p x - A special emphasis was 
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FIG. 6: Confidence regions in the u — a plane according to 
the sample of angular size data by Gurvits et al. [i^ ] and 
fixed Qm — 0.263 as shown in panel. The confidence levels of 
the contours are indicated. The point "x" marks the best fit 
values, to — —1.03 and a = 0.90. 



FIG. 7: Confidence regions in the £Im — ot plane according to 
the sample of angular size data by Gurvits et al. !42l| . For a 
phantom cosmology with uj = —1.023, the confidence levels of 
the contours are indicated. As in Fig. 6, the "x" also points 
to the best fit values shown in the panel. 



given to the case of phantom cosmology (to < —1) when 
the dominant energy condition is violated. The effects 
of the local clustered distribution of dark matter have 
been described by the "smoothness" phenomenological 
parameter a(z), and a simple argument for its functional 
redshift dependence was given in the Appendix A (see 
also Figure 1). 

The influence of the dark energy component was quan- 
tified by considering the angular diameters for sample 
of milliarcsecond radio sources (Fig. 5) as described by 
Gurvits et al. [52] ■ By marginalizing over the character- 
istic angular size /, fixing f2jvf = 0.263, and assuming a 
Gaussian prior on the EOS parameter, i.e., u> = — 1±0.3, 
the best fit model occurs at u> = —1.03 and a — 0.9. This 
phantom model coincides with the central value recently 
determined by the Supernova Legacy Survey (Astier et. 
al. Q). On the other hand, fixing ui = —1.023 and as- 
suming a Gaussian prior for CIm , that is, 0,m = 0.3±0.1, 
we obtained the best fit values (f^f = 0.29, a — 0.9). 

Finally, in order to improve the present results, a sta- 
tistical study is necessary for determining the intrinsic 
length of the compact radio sources. Further, unlike to 
what happens with SNe data [4l|, the angular diameter 
sample of compact radio sources of Gurvits et al. [42j 
does not provide useful constraints on the a parameter 
(see Figs. 6 e 7). Naturally, these results reinforce the 
interest for measurements of angular size from compact 
radio sources at intermediary and high redshifts in order 
to constrain the a parameter with basis on the ZKDR 



distance. 



APPENDIX A: ON THE REDSHIFT 
DEPENDENCE OF a(z) 

In this Appendix we discuss the functional redshift de- 
pendence of the smoothness parameter, ev(z), adopted in 
this work. By definition 



ot(z) = 1 — 



Pcl(z) 

Pm(z) ' 



(Al) 



where p c \ denotes the clumped fraction of the total mat- 
ter density, pm, present in the considered FRW type Uni- 
verse. This means that the ratio between the homoge- 
neous (ph) and the clumped fraction can be written as 
Ph/ Pel = a ( z )/[l ~ a ( z )]- How this ratio depends on the 
redshift? In this concern, we first remember that a(z) 
lies on the interval [0,1]. Secondly, in virtue of the struc- 
ture formation process, one expects that the degree of 
homogeneity must increase for higher redshifts, or equiv- 
alently, the clumped fraction should be asymptotically 
vanishing at early times, say, for z > 100. This means 
that a(z) — > 1 at high z. On the other hand, a must 
be zero for a completely clustered matter which is dis- 
proved at low redshifts by the data of galaxy clusters [H . 
It thus follows that at present [z = 0), the related frac- 
tion assume an intermediate value, say, (3 Q . In addition, 
it is also natural to suppose that the redshift dependence 
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of the total matter density, pm , must play an important 
role in the evolution of their fractions. In this way, for 
the sake of generality, we will assume a power law 



Pk_ a{z) _ ,Pm, 7 

— i i \ — An I ■ 
Pd 1- a(z) p 



(A2) 



where f3 a = {ph/ Pci)z=o and 7 are dimensionless num- 
bers. Finally, inserting pm{z), and solving for a(z) we 
obtain: 



a(z) = 



/3o(l + z) 37 
1 + /3 (1 + z) 3 t' 



(A3) 



which is the expression adopted in this work (see Eq. 

(hod). 

As one may check, for positive values of 7, the smooth- 
ness function (A3) has all the physically desirable proper- 
ties above discussed. In particular, the limit for high val- 
ues of z does not depend on the values of (3 and 7 (both 



of the order of unity). Note also that if the clumped 
and homogeneous portions are contributing equally at 
present (/3 = 1), we see that a(z = 0) = 1/2 regardless 
of the value of 7. Figure 1 display the general behav- 
ior of a(z) with the redshift for different choices of j3 
and 7. The above functional dependence should be com- 
pared with the other ones discussed in the literature (see, 
for instance, [13; liSl [22( and Refs. therein). One of the 
most interesting features of (A3) is that its validity is not 
restricted to a given redshift interval. 
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